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Abstract 

Recently there has been theoretical and experimental interest in Bloch-Siegert shifts in an intense 
photon field. A perturbative treatment becomes difficult in this multiphoton regime. We present a 
unitary transform and rotated model, which allows us to get accurate results away from the level 
anticrossings. A simple variational energy estimate leads to a new expression for the dressed two- 
level system energy which is accurate, and useful over a wide range of the dimensionless coupling 
constant. 

PACS numbers: 32.60.+i,32.80.Bx,32.80.Rm,32.80.Wr 



1 



I. INTRODUCTION 



The basic problem of a two-level system coupled to an oscillatory perturbation arose in 
the early on in quantum mechanics in association with studies of spin dynamics in a magnetic 
fieldfii^. The closely related problem of a two-level coupled to a simple harmonic oscillator 
has been of interest for more than 30 years, at least since the work of Cohen- Tannoudji et 
al^. The Hamiltonian for this problem can be written as 



This model has been of interest recently in problems in which atoms interact with strong 
electromagnetic fields. We have become interested in this problem since it can exhibit 
coherent energy coupling between systems with characteristic energies that are very different. 
The multiphoton regime in which AE ^ hujo in this regard is of interest to us. Additionally, 
our focus has been on the large n limit. 

It is well known that the dipolar interaction between photons and atoms causes the 
frequency of the atom to get shifted. The associated shift has been termed the Bloch- 
Siegert shift in the literatureii^. A general analytic expression for the {2k + l)th resonance 
has been known for some time^. Measurements of this shift have been reported recently up 
to k = 11-. Various theoretical methods have been used to understand these results^i^. 

We have recently found a unitary transformation that separates the problem into three 
parts, each of which is relatively complicated. In the multiphoton regime we view the 
largest of these as an "unperturbed" dressed Hamiltonian, which we have found give energy 
eigenvalues that are close to those of the original problem away from the level anticrossings. 
In the present work, our attention is drawn to this aspect of the problem. Using this 
approach, we are able to develop a reasonably good estimate of the Bloch-Siegert shift using 
a simple variational approximation, one which we can compare with results obtained using 
other methods. 

II. ROTATED HAMILTONIAN 

As mentioned above, we have found it helpful to work in a rotated frame. Instead of 
considering the Hamiltonian in Equation ([1]) , we consider the unitary equivalent Hamiltonian 



(1) 



where we can write the spin operators Sj in terms of the Pauli matrices as 
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where A is 
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The operator A depends explicitly on a + a) . 

This rotation can be carried out explicitly, leading to a complicated expression. It is 
useful to write the rotated Hamiltonian in the form 



H' = Ho + V + W 
where we have defined Hq, V, and W according to 
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Even though this rotated Hamiltonian is much more complicated than the one we started 
out with, we have found that it has useful properties that are of interest. For example, from 
numerical calculations and also from the WKB approximation, we have found that the "un- 
perturbed" part of the Hamiltonian is a good approximation to the initial Hamiltonian 
H away from level anticrossings. The "perturbation" V contains the part of the Hamiltonian 
that is most responsible for level splitting at the anti-crossing (as we will discuss elsewhere). 
Finally, the remaining term W gives rise to a minor correction to Hq. 



III. BLOCH-SIEGERT SHIFT 

To illustrate the usefulness of this rotation, we focus on one issue in particular: the devel- 
opment of an interesting approximation for the energy eigenvalues. For this, we need focus 
only on Hq. Consider the time- independent Schrodinger equation for the "unperturbed" 
part of the rotated problem 

^^V' = ^oV' (7) 
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Because Hq is somewhat complicated, it is not straightforward to develop exact analytic 
solutions. We have used numerical solutions to develop some understanding of the eigen- 
functions and associated eigenvalues. From these numerical studies, it has become clear that 
there exists a useful analytic approximation. 
We adopt an approximate solution of the form 

ipt = \n)\s,m) (8) 

where \n) is a pure SHO eigenstate, and \s,m) is a spin function. The variational estimate 
for the energy eigenvalue in this approximation is 

This approximate energy can be evaluated explicitly to give 



(9) 
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This result is interesting in that it leads directly to a simple interpretation of the system 
in terms of a dressed two-level system and unperturbed oscillator, which is reasonably good 
away from the level anticrossings. Moreover, the dressed two-level system energy is involved 
in the calculation of the Bloch-Siegert shift in the literature. 
Let us write the variational energy as 



where we define g according to 



Et = AE{g)m + hwQn 



AE 
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and AE{g) as 
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This approximation is very good, as can be seen in Figure [T] where the approximation is 
plotted along with the exact numerical results. The Bloch-Siegert resonance condition can 
be written in terms of AE{g) as 

AE{g) = {2k + l)huJo (13) 

It is possible to develop a simpler approximation for this approximate dressed energy 
based on the WKB approximation. We may write 
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FIG. 1: Plot of AE{g) computed from direct numerical solutions for AE = llhujQ and n = 10^, 
and also the approximate AE{g) based on the approximation of Equation (jl2p . The two curves 
are so close that they cannot be distinguished. 
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AE{g) = -- / J '-^dy (14) 



with e = 2n + 1. Results from this approximation can also not be distinguished from the 
exact results if plotted as done previously in Figure 1. 



IV. COMPARISON WITH OTHER RESULTS 



As the results available in the literature of interest here are given for small g, we require 
a power series expansion of Equation (|T2l) : 
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To compare to the exact expression obtained by Ahmad and Bullough^, we write the 
resonance condition using their notation as 
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We translate this into our notation to obtain 
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Keeping in mind the resonance condition [Equation f[T^ ]. this is consistent with 
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A^(^) = l + a{k)g'^ + b{k)g^ + ■■■ 
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where 
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We can now see that in the large k hmit, Equation (ITSll reduces to our result. It should also 
be noted that the deviations for finite k are 0(l//c^). 

Another approximate multiphoton result is given by Ostrovsky and Horsdal-Pedersen^. 
This result (written in their notation) for the Bloch-Siegert shift is 



This result is equivalent to our result [Equation f|T5l) ] to fourth order in g. 
V. CONCLUSION 

The unitary equivalent Hamiltonian provides a simple way to understand the Bloch- 
Siegert shift. We have found in our work that the unperturbed part of the rotated Hamil- 
tonian [Hq) produces a reasonably accurate estimate of the dressed two-level system energy 
in the large n limit away from the level anticrossings. In the discussion presented here, 
we have shown that a simple variational estimate for the energy associated with Hq in the 
rotated frame produces good results away from the level anticrossings over a wide range of 
the dimensionless coupling strength g. The approximation is found to be in agreement with 
previous work in the limit that g is small. 
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